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Abstract 

We show that gaussian quantum fluctuations, even if inflnitesimal, are suf- 

flcient to destroy the superfluidity of a disordered boson system in ID and 

2D. The critical disorder is thus flnite no matter how smaU the repulsion is 

between particles. Within the gaussian approximation, we study the nature 

of the elementary excitations, including their density of states and mobility 

edge transition. We give the gaussian exponent ry at criticality in ID and 

show that its ratio to r] of the pure system is universal. 
PACS numbers: 67.40.Yv, 74.20.Mn, 05.70.Jk, 75.10.Nr 
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In the presence of disorder, (repulsive) interacting bosons can undergo a transition from 
the superfluid (SF) phase into an insulating Bose-glass (BG) phase IQ - [if- This transition 
is intrinsically quantum in nature in that no amount of disorder will destroy the superfluidity 
without invoking the non-commutativity of density p and phase 0. Hence, the usual saddle 
point or Hartree solution is always long-range ordered, and corresponds to a non-uniform 
condensate. Given that, it is of interest to investigate what 'minimal' quantum effects are 
necessary to give a transition. In terms of going beyond the saddle point approximation, 
these effects can be characterized as gaussian, non-linear, topological (as in vortices) etc. 
In effect, one is asking 'what drives the transition', even if the true universality class of 
the transition may require quantum fluctuations beyond the 'minimal' ones. To clarify this 
perspective, consider the 2D classical XY model as an analogy. There, the true long range 
order (LRO) is destroyed at any finite temperature by spin waves, even though (bound) 
vortices do renormalize the exponent rj (universality class). That is, spin wave alone can 
explain why the low temperature phase has algebraically decaying correlations. On the other 
hand, vortices must be invoked to explain the Kosterlitz-Thouless transition [Q. 

In this article we show that gaussian fluctuations, even if infinitesimal, are sufficient to 
destroy superfluidity in ID and 2D at finite disorder. The model we use is the hard-core 
boson model with on-site disorder, which is equivalent to the spin-1/2 XY magnet with 
a transverse random field Written in a rotated frame for later convenience, the 

Hamiltonian is 0: 

^ = E (^r^l + ^r^i ) - E h,sy , (1) 

<i,3> 3 

where the random field {hj} is given by independent gaussian distribution function P{hj) 
with width h. We pick this model because while it contains the features believed to be 
essential for the SF-BG transition, it has a simple classical solution and a 'built-in' parameter 
to systematically investigate quantum fluctuations. It is also of interest as a disordered 
quantum spin system 0-10. Our study reported in this paper will be focused on the 
T = case. 
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The off-diagonal LRO of the boson system is related to the magnetic LRO in the x — z 
plane. The classical solution for this model corresponds to treating the spins as classical 
vectors. In terms of the bosons, one is dealing with variational wavefunction of the form 
Y[j{uj + Vjbj')\0), which when projected into states of definite N are Jastrow wavefunctions 
given by Gutzwiller projection of condensate wavefunction -^b'j)^\0) . In Ref. it 
was shown in the spin model that provided there is no gap in the spectrum, the classical 
ground state is always ordered. However, it was also seen that the ground state is not long 
range ordered with strong disorder when the quantum (specifically S = 1/2) nature of the 
spin operators is taking into account. Following the motivation discussed in the previous 
paragraph, it is thus of interest to investigate if the destruction of LRO can be achieved by 
gaussian quantal effects. This can be studied by means of a spin wave analysis ||^. 

Within such an approach, the first question is what would be the signature of destruction 



of the LRO Since the spin-wave analysis is an expansion about the ordered state, this 
destruction is indicated by an instability. Possible scenarios are 1) a diverging fluctuation in 
the order parameter, 2) negative excitation energies, or 3) complex excitation energies (e.g., 
Bogoliubov's solution to bosons with attractive interactions). In the pure case, scenario 1) 
is observed in ID. Exact solution for S = 1/2 |TT|] and general understanding of ID spin 



systems indicates that this diverging fluctuations destabilize the LRO and the ground state 
has algebraically decaying correlation functions. 

We now derive the spin- wave Hamiltonian 0. First we generalize Hamiltonian (^ to 
arbitrary spin S by rescaling J J / S"^ and hj hj/S. In the inflnite S limit, the spins 
behave classically. Taking the z-axis as the ordering axis, the spin on site j lies on the y — z 
plane at angle 9j from the z-axis, with {9j} given self-consistently by 

sinOjJ ^ cos9j' = hj cos 9 j , (2) 

<j'> 

where < j' > indicates nearest neighbors of the site j. The statement that LRO persists to 
all order is revealed by the solution to (|^) having all cos 9j 7^ no matter what value of h is. 
A local rotation about the a;-axis is performed so that the spin points along the new 2;-axis. 



The usual Holstein-Primakoff transformation |T2[ of the spins into boson operators can now 
be defined in the rotated frame. To order 1/S, one arrives at a quadratic Hamiltonian for 
the bosons P: 

H = - J cos6'i cos6'j - Y sin6'j - y^; XI {jijO^Uj + Kija^a^ + H.c)j + 0{-^) , 

<i,3> j <i,3> 

(3) 

where Jij = J(l + sin^jsin^j) + -^^Sij and Kij = J(l — sin 6^, sin ^j), which describes 
gaussian fluctuations of strength 1/5* about the classical ground state. In Ref. 0, (0) is 
studied perturbatively for weak disorder. In this paper, we will diagonalize numerically 
on finite-sized lattices, and will not limit ourselves to weak disorder. This will enable us to 
study the destruction of LRO. (|) is formally diagonalized by a Bogoliubov transformation 

«i = II(Wi«7a + ^^ia7i) , (4) 
a 

where a is the eigenstate index. We have taken the u's and f 's to be real. The 7's are boson 
operators if 

X](MjQ,MjQ,' — VjaVja') = ^aa' ; (5) 

i 

and we seek the solution 

n3W = E' + Y^a7l7a , (6) 

a 

which implies the Bogoliubov equations for m's and f 's. 
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Y {Jjj'Uj'a + KjjiVj'a) , 



^'"ja = Y i^jf^j'a + Jjj'Vj'oc) , (7) 

<j'> 

to be 'normalized' by the condition (^. For sites, this is a 2N x 2N matrix equations with 
2N eigenstates. Note that for a given solution with eigenvalue u, there is the complimentary 
solution u ^ V, with eigenvalue —u. However, only one of these can be consistent with @, 



and the other is unphysical, leaving us with physical solutions. The Goldstone mode, 
corresponding to uniform spin rotation about the y-axis in (|lj), is given by Ui = Vi cc cos 6i. 

We investigate LRO instability in ID and 2D. Calculations in ID are done on lattices 
of size 50 - 120, averaging over 500 configurations for each value of A = J/h, and in 2D 
on6x6tollxll lattices averaging over 200 configurations. Instability criteria 2) and 
3) are not observed, leaving 1), a diverging fluctuation in the order parameter as the sole 
possibility. Within the spin wave approximation as formulated, the relevant quantity is 

= ^ E cose,S{S]) = ^ E E ^ose.vl = J du;N{uj)v\u;) , (8) 

j j aj^O 

where N{lj) = jj: Y,a ~ ^a) is the density of states (DOS). 

As remarked earlier, in ID 6m diverges as — oo even without disorder, so it seems 
criteria 1) is inapplicable. However, more precisely, 6m oc InA^, and we view this as an 
indication for an algebraic LRO (replacing the true LRO, see later), hence the ground state 
is still a superfiuid. Thus, we argue that the transition is marked by 6m diverging faster 
than InA^. This is in fact seen in our calculation, and is shown in Fig. 1, with the critical 
value of A = Ac ~ 0.6 in the present model. The transition occurs thus at finite disorder. 
Since the spin wave approximation is correct in the large S limit, there is a discontinuity 
between S ^ oo and S = oo. In 2D, 6m/ m is finite in the pure case as A^ — » oo, which we 



take to mean the LRO is stable, and is consistent with the exact result for S* = 1/2 |jT^. Fig. 
2 shows 6m/ m vs. InA^ for different values, and we see that there is a transition between 
A = 0.1 and A = 0.08. 

Thus, already in the gaussian approximation, in contrast to the classical case, there is a 
transition from an algebraic in ID and a true long range ordered in 2D superfiuid phase to 
a disordered phase. Since in a gaussian theory, the ground state is just the classical state 
modified by the zero point motion of the the excitations, it is of interest to ask whether the 
transition is due to a change in the DOS {N{uj)) or the nature of the excitations {v'^{uj)) 
or both. In the pure case, f^(u;) oc ^ for small tu, while N{ijj) oc uj'^~^ for small uo. For the 
infinitely strong disorder (J = 0) case, the excitations are single spin flips, with excitation 



energies \hj\. Hence N{u!) is simply given by the distribution of hj, and is finite at low 
energies. It seems reasonable to expect therefore Nq — N{(jj 0) is finite in ID for all 
A, and the transition must be due to ^^(a;) diverging faster than l/cu. This picture is 
confirmed by our numerical calculations and Fig. 3 is shown for DOS in ID. While there 
is some ambiguity in deciding Nq for infinite system from finite-size calculation, we have 
checked to see that the scaling of A^o with N is in fact consistent with a non-zero DOS at 
zero energy. For A < Ac, Sm oc N^, with 9 — 0{A.). We find that this exponent is in 
agreement with the exponent of v'^{u) oc with 5—1 + 9, again indicating A^o finite. In 
2D, the ordered phase should be characterized by N{u!) oc lu and the disordered phase by 
A^o finite. Our results are consistent with this. For u; > 0.1, N{u;) is linear in a;, with the 
slope increasing with decreasing A. For A < 0.08, A^o is finite. Unfortunately, we cannot 
say for certain whether the DOS transition exactly occurs at the order parameter transition 
due to the inability of pinpointing Ac (The popular finite-size scaling method for locating 
critical point is not applicable here since there is no scale invariance). We hope to clear up 
this point in a later publication. 

The low-energy excitations calculated here are particularly significant in the SF phase, 
as they are approximations to the collective modes (phonons). These excitations can be 
extended or localized. It is of interest to ask if their localization transition is related to the 
'locahzation' of the ground state. It is also of interest by itself as an Anderson localization 
problem of the eigenstates of (3). Since the zero-mode corresponds to uniform phase rotation, 
it must be extended. One thus expects that possibly for a given A, a transition from 
extended to localized states with increasing energy at a mobility edge energy Ec, and perhaps 
as A — > Ac+. The quantity we calculate as a measure of localization is the inverse 
participation ratio p, which we assume for localized states scale as the inverse localization 
length and is zero for extended states in infinite systems. Ec is the energy where p first 
vanishes. However, for finite size, p scales as the greater of L"^. Hence, at low energy, 
where ^ > L, p{E) is constant, and only for E > Ec{L), where p{Ec{L)) = L~^, does p{E) 
gives the behavior of an infinite system. One way to obtain Ec is by extrapolating the part 
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of the curve to where p = 0. An improved method is to extrapolate Ec{L) to L — oo |T5 . 
Eq. d^) guarantees that Uja and Vja are either both extended or locahzed, so it suffices to 
calculate p{E) for uja- In Fig. 4 we show p{E) for various L's in ID for A = 1.5 (SF side) 
and A = 0.5 (insulating side), which for the former seems to show clearly a finite Ec. For 
the latter, we ascertain E^ to be very small, probably zero (our extrapolation actually gives 
an unphysical small negative value). This seems to support the idea that the localization 
of the ground state and the excitations occur simultaneously, a feature of Giamarchi and 
Schulz's theory in ID [|T^] . 



Upon refiection, however, we have serious doubts. In the perturbative (in disorder) 
calculation of Ref. ||^, the phonon mean free path is found to diverge as E~^'^'^^\ Common 
wisdom has it that in ID, the localization length and the mean free path are essentially 
identical, since any scattering is backscattering. Hence, one expects the relatively slow 
C,~^{E) oc E"^ for small E, crossing over to a more rapid E dependence at a higher cross-over 



energy E^. This is in fact known to be the case for classical vibrational modes in ID ||17|| , 
and all states except the uniform translation mode are localized. The problem is then that 
the divergence of the localization length ^ is not a true critical phenomena with a critical 
exponent. In our calculation, if the size L < C,{Ex), then we cannot probe the weak E 
dependent regime, and we will mistake E^ for the true mobility edge. For comparison, we 
look at a finite system of random masses connected by springs and find p{E) curves similar 
to Fig. 4. While our results do not constitute evidence for all eigenstates of (^ to be 
localized in ID, we believe this is in fact the case, and what Fig. 4 shows is E^ decreasing as 
the disorder is increased, vanishing at or close to the superfiuid-insulator transition. Since 
even in 2D, it is believed that all classical waves are localized it seems probable that 
all phonon excitations are localized too. This fact of localization of all the excitations, if it 
is true, may invalidate the usual effective field theory based on the action of propagating 
phase modes, which is crucial for the scaling theory of Ref. 0. Indeed, certain predictions 
of the scaling theory has been questioned by recent quantum Monte Carlo simulations 
on 2D hard-core dirty bosons. 



In our model, the classical state is the Gutzwiller state, and 1/S serves as an expansion 
parameter for quantum fluctuations. We find a transition at finite disorder when only 
gaussian fluctuations are kept. Since the physics is sufficiently general, we believe this 
conclusion would hold true if one consider soft-core (e.g., Hubbard U) models and use fi as 
the quantum expansion parameter. What about a phase diagram of disorder vs. Ul Since 
bosons will condense into the lowest energy (localized) state for [7 = 0, the critical disorder = 
0. How about f/ — » 0? The problem is that increasing U both affects the classical condensate 
and enhances quantum fluctuations. In fact, the Hartree solution becomes extended with 
infinitesimal U . Thus, the f/ ^ limit should not be qualitatively different from the limit 
of 1/S 0, and the critical disorder is again finite. This conclusion is in agreement with 
numerical works performed on the disordered boson Hubbard model 

In ID, Ref. [1^ predicted that the renormalized critical exponent rj is universal and equal 



to 1/3 at the SF-BG transition based on a perturbative renormalization group calculation. 
One might ask what the value of rj is in the gaussian theory. Rigorously, the spin wave 
theory as formulated cannot produce a power-decaying correlation function (it is similar 
to using (Vtt)^ as the action in the classical non-linear a model ). However, r] and 7, the 
coefficient of In in 5m/ are proportional in the pure case. Assuming the relation holds 
even with disorder, it implies 

(9) 



Vpure ''ipure 

From the slope of the A = 0.6 curve in Fig. 1, we estimate rjc/rjpure ~ 1-4. This ratio is 
universal, while rjc is not. ?7c ^ as l/S* — in our model, or as f/ — in Hubbard type 
models. 
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FIGURES 

FIG. 1. Fluctuation corrections to the order parameter 6m /m is plotted against IniV. The 
divergence becomes faster than In as A exceeds a critical value Ac ~ 0.6. The solid lines are 
obtained through a linear fit and the dashed line is a guide to the eyes. 

FIG. 2. The same as Fig. 1, plotted for 2D systems. Unlike it in ID, 5m is finite for weak 
disorder and diverges for A > Ac, which is between 0.1 and 0.08. (6) is the same plot as it in (a), 
but is presented on a different scale, which shows clearly that 5m is bounded as AT ^ oo for small 
A. 

FIG. 3. DOS on the insulating side (A = 0.5) in ID. Here N = 100. 

FIG. 4. Participation ratio p{E) in ID is plotted for several values of N. The pseudo-mobility 
edge Ec is obtained from these plots through the procedure described in the main text. In (o), 
A = 1.5, the system is in the superfluid phase and Ec is finite. Ec is about zero in the insulating 
phase as it shows in (6) for A = 0.5. 
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